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Abstract 

Beisert et al. have identified an integrable SU(2, 2) quantum spin chain which gives 
the one-loop anomalous dimensions of certain operators in large N c QCD. We derive a 
set of nonlinear integral equations (NLIEs) for this model, and compute the scattering 
matrix of the various (in particular, magnon) excitations. 
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1 Introduction 



The search for integrability in QCD has a long history (see e.g. p]-[6] and references therein). 
A remarkable recent development is the discovery [7J that the one-loop mixing matrix Q for 
the chiral gauge-invariant operators 

Pi ( x ) ■ ■ ■ fa L /3 L {%) (1-1) 

in the limit N c — > oo is given by the integrable spin-1 antiferromagnetic XXX Hamiltonian 

mm, 



L r 



Z=l 



7 1 - - 1 - - 



1.2) 



Here f a/ 3 are the selfdual components of the Yang-Mills field strength F^ u = d^A v — d u A^ — 
igYMiA^ , A u ] (where the gauge fields A^x) are N c x N c Hermitian matrices), which together 
with the anti-selfdual components are defined by 

where = ia 2 (a tl a v -a v a tJ )/4, = -i(a li a v -a v a tl )a 2 /4 : and cr M = (1 ,a), a M = (1 , -a). 
Moreover, a s = g\ M jAm^ a s N c is the 't Hooft coupling [1] which is assumed to be small, and 
5* are spin-1 generators of SU{2). Indeed, since f a p has three independent components 

/+ = /!!, /0 = ^(/l2 + /2l), /- = /22, (1-4) 

the operators fll.ip can be identified with the Hilbert space of a periodic spin-1 quantum spin 
chain of length L. The eigenvectors and eigenvalues of T, i.e. the linear combinations of the 
operators f ll.ip which are multiplicatively renormalizable and their anomalous dimensions, 
respectively, can therefore be obtained using the Bethe Ansatz [TTI [T2]. In particular, the 
anomalous dimensions are given by 

7 = ^I^-V7^1 • n-> 




2vr 

where {l\ , . . . , /m ; } are roots of the Bethe Ansatz equations (BAEs) [§ 



1 Given a set of operators O m (x), the mixing matrix is defined by Y = Z^ 1 ■ dZ/dlnA, where Z is 
the renormalization factor which makes correlation functions of 0~. n (x) = Z a ^O n (x) finite, and A is the 
ultraviolet cutoff. See also [8]. 

2 There is an additional (zero-momentum) equation due to the cyclicity of the trace in the operators. 
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This result was generalized in [13J to gauge- invariant operators with derivatives 



tr {D m > /) . . . (D m <- f) , (1.7) 

where 

D m f = £> ai6l . . . D olm6lm f Pl + symmetrized (1.8) 

(complete symmetrization in the undotted and dotted indices, respectively), and = 
a® a D a a is the usual Yang-Mills covariant derivative. Namely, the one-loop mixing matrix for 
the operators (11.71) is given by an integrable 50(4,2) = SU {2,2) (non-compact!) quantum 
spin chain Hamiltonian with spins in the representation with Dynkin labels [2,-3,0]. The 
anomalous dimensions are given by 

7 = ^|^->; 7 ^T + >;zrT : 7w7) , (L9) 




where the BAEs are now given by 
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I j + M -TT lj - k + \ TT lj ~U k - 1/2 

h-i) Ll h - l k-i ^ lj -u k + i/2 ' 

Uj - 3l/2 \ = yr Uj-Uk + l yr Uj - l k - l/2 yr Uj - T k - l/2 

uj + 3i/2 J Al u . - Uk -i Uj -l k + i/2 1 = 1 U j - r k + i/2 ' 

A/ r . ■ M u . ,~ 

TT r j ~ ^fc + » TT r j - Mfc - V 2 

As noted by Beisert et al., a w-root corresponds to adding a covariant derivative D^; and 
an /-root and an r-root flip a left-Lorentz-spin 1 — >• 2 and a right-spin 1 — > 2, respectively. 
The scaling dimensions and SU{2) l x SU {2)r quantum numbers are given by 

D = 2L + M U , S 1 = L + ^M U -M l) S 2 = ^M u - M r , (1.11) 

respectively. 

As noted in [13], the BAEs (11.101) can be obtained from those of the "beast" form of 
Af = 4 SYM [TJ] by truncating the supergroup SU{2,2\4) down to the Bosonic subgroup 
577(2,2). § Much attention has been focused on the S matrix of Af = 4 SYM and of the 
corresponding string theory (see e.g. (16] ). 



3 For some early references on integrable gl{n\m) spin chains, see e.g. |15j 
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For the pure spin-1 problem (11.51) . (11.61) . the ground state for large L is described by 
a "sea" of approximate "2-strings" of Z-roots [UJ [12] (in contrast to the case of the spin- 
1/2 antiferromagnetic XXX chain, for which the ground state is described by a sea of real 
roots). The excitations consist of "spinons" (roughly speaking, "holes" in the sea) which 
carry RSOS [17] quantum numbers. The spinon-spinon 5* matrix was found by indirect 
methods in [TBI EE], correcting the result obtained in [TTJ using the string hypothesis. A 
nonlinear integral equation (NLIE) [201 EI] has been obtained for this model [22]-[2I], which 
does not rely on the string hypothesis and provides a more direct way to compute the S 
matrix [25]. The NLIE of the SU(2) sector of M = 4 SYM has been studied in [26]. 




For the general case (jl.9p . (11.101) . the ground state is still a sea of approximate 2-strings 
of /-roots, since the w-roots contribute positively to the energy (and the r-roots do not 
contribute at all). Hence, there are again spinon excitations corresponding to holes in the 
sea. However, there are now also "magnon" excitations, corresponding to u- roots [T3] . 

Our main objective here is to further investigate these magnon excitations, and in par- 
ticular, to compute the magnon-magnon 5* matrix. Owing to the nontrivial nature of the 
ground state, this S matrix (like the spinon-spinon S matrix) must be computed with care: 
using the string hypothesis as in [UJ gives an incorrect result. To this end, we first derive in 
Sec. |2] a set of NLIEs for the model. Although we do not invoke the string hypothesis, we 
do make a certain analyticity assumption in order to describe the w-roots. For simplicity, 
we restrict to real -u-roots, and we do not consider r-roots. We then use these NLIEs to 
determine the energy and momentum of the excitations (Sec. [3]), and their S matrices (Sec. 
HJ) . We end in Sec. [5] with a brief discussion of our results. 

2 Nonlinear integral equations 

We restrict our attention to the case without r-roots (M r = 0), for which the BAEs (11.101) 
reduce to 






- TT j ~ Uk ~ % J— 



(2.1) 



Uj — 3i/2\ L 
Uj + 3i/2J 



1V1 U 

n 




Uj — Uk + i -i-i- Uj — Ik — i/2 



Uj -u k -i J-jL Uj -l k + i/2 ' 



(2.2) 



We now proceed in turn to recast these two sets of BAEs in the form of NLIEs. 
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2.1 The first set of BAEs ( 12.11) and an auxiliary inhomogeneous 
mixed spin chain 



An important hint on how to analyze the first set of BAEs (12. ip comes from rewriting it in 
the obviously equivalent form 

lj + i\ L yr lj - u k + i/2 _ j± lj - l k + i 

11/ ._,„ _,-/9 "11/ -7, -J" [Z - 6) 



h-i) t=i l i~ Uk ~ i / 2 



k=l 3 



k 



We recognize these as the BAEs for an inhomogeneous "mixed" spin chain which has two 
types of spins: spin-1 and spin- 1/2, with L of the former and M u of the latter. (See, e.g., 
[27].) Moreover, the latter have associated "inhomogeneities" iuk, k — 1, . . . , M u . 

We therefore consider an auxiliary integrable inhomogeneous mixed quantum spin chain, 
where the number of spin-1 and spin- 1/2 "quantum" spaces are given respectively by L and 
M u ; and with spectral parameter inhomogeneities iuk only for the spin- 1/2 spins. This chain 
has two relevant transfer matrices 7\ (x) , T 2 (x) , corresponding to "auxiliary" spaces which 
are spin-1 /2 (2-dimensional) and spin-1 (3-dimensional), respectively. 

We find by standard methods that the eigenvalues of these transfer matrices (which we 
denote by the same notation) are given by0 

Ux) = i ;( x - i /2) ( p(x- i )^^ + ^(x + l /2) ( p(x + i )^^-, (2.4) 
T 2 {x) = 4>{x)^{x - i)<f>(x - i/2)<j>{x - 3i/2) Q ( X + 3 ' /2) 



+ i^{x) 2 (j)(x- i/2)(j)(x + i/2) 



Q{x-i/2) 

Q{x + 3i/2)Q{x-3i/2) 
Q{x + i/2)Q(x-i/2) 

Q(x - 3i/2) 



+ iPtx^x + i^x + i^^x + Zi^) . 

Q[x + 1/2) 

■= Ai(x) + \ 2 {x) + X 3 (x) , (2.5) 



where 



M u Mi 

<P(x) = x L , ^(x) = H(x- Uj ), Q(x) = Y[(x-l j ). (2.6) 

i=i j=i 

Indeed, the BAEs obtained by demanding that Ti(x) be analytic at x — lj (zeros of Q(x)) 
coincide with ( 12. lft . 



Note that in place of the standard spectral parameter u, we introduce u = ix. 
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Evidently T 2 (x) has the common factor ip(x), which has "trivial" zeros. We therefore 
introduce the renormalized T 2 , 

T 2 (x) = ^{x)t!?\x). 

We note that 

S 1 = S z 1 =L+^-M l , (2.7) 
and we recall that the "energy" is given by ( 11.91) . 



2.1.1 Physical degrees of freedom 

For simplicity, we restrict Uk to be real, and M u = 1,2. For now, we also assume that are 
given by hand, with (in the case M u — 2) u\ — —u 2 . We shall discuss how they should be 
determined later in Sec. 12.21 

Numerical studies for small values of L suggest that: 

• For M u = 1, the lowest energy state in the Sf = 1/2 sector is characterized by a single 
zero (^q.) of T\(x), and a single zero (6^) of T 2 \x). Both of these zeros lie on the real 
axis. 

• For M u = 2, the lowest energy state is in the Sf = sector. In the "physical strip" 
(—1/2 < Smi < 1/2), T\(x) and T 2 \x) are free from zeros. 

• For M u = 2, the second- lowest energy state is in the = 1 sector. It is characterized 
by two zeros ($ a ) of Ti(x) and two zeros (Oh) of T 2 \x). These zeros lie on the real 
axis. 

These observations suggest that three sets of real parameters are needed to describe the 
physical degrees of freedom: Uj ,$ a ,8h- The first and third parameters correspond to magnon 
and spinon rapidities, respectively. The second parameter, which seems to correspond to 
excitation of the RSOS degree of freedom, is not discussed in [13] . 

2.1.2 The auxiliary functions and algebraic relations among them 

As in previous studies [221 121], we introduce a pair of auxiliary functions 
U f\ Ai(x) + A 2 (x) - X 2 (x) + X 3 (x) 

bi(x) := -— ismx>0, oAx) := r~T~\ ismx<0, (2.8) 

X 3 (x) XAx) 
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where Xi(x) are defined in (12. 5p . They are free from zeros and poles near the real axis. This 
will be apparent from the following representations, 

l[) ^<f> 3/2 <f> l/2 Q(x-3i/2) ll[X h U) > 

= . / 1/2 , g^§ r l(g+ i/2). (2.9) 

^-10-3/20-1/2 Q{ x + 3V 2 ) 

We have introduced here the abbreviated notation <p a := <p(x + ia), and similarly for ip, which 
we shall use throughout this part of the paper. 

At this stage, there seems to be no reason why the two auxiliary functions should be 
introduced in the corresponding half planes. This will become clear at a later stage. 

The upper-case functions are also introduced: Q3i(x) = 1 + b±(x), %$i(x) = 1 + bi(x), 
and the following relations are also useful: 

7f(x) = Q^y^ BiW (2.10) 

- ^-^Wm %lix) - (211) 

Apparently *Bi(x) vanishes at x — 6h, but it remains nonzero at x — Uj. 
We now define the most important functions, 

b(x) = bi(x + ie) , 35 (x) = »i(ar + ie) 3mi>0, (2.12) 

b(x) = bi(x - ie) , »(z) = » x (a: - ie) 3mi<0. (2.13) 

Here e denotes a positive quantity which is slightly larger than the deviation of the 2-strings 
from their "perfect" positions. Therefore 23 would possess zeros (due to the factor Q(x + i/2) 
in ( 12.10(1 ) slightly below the real axis if it were defined in the whole complex plane. The 
function 03 is, however, defined only in the upper half plane (including the real axis). 

Another auxiliary function originates from the so-called fusion formula that relates the 
two transfer matrices, 

T 1 {x-i/2)T 1 {x + i/2) = ^^-103/20-3/2 + ^oT 2 (r) (a;), (2.14) 

which can be verified using (12.41) and ( 12.5(1 . For later convenience, we renormalize T\(x) = 
Yla=i tanh ~{x — i? a ) T± \x), and rewrite the above in the form 

T[ r) (x - i/2) T[ r) {x + i/2) = ^1^-103/20-3/2 + ^oT 2 (r) (x) 

= ^ 1 ^_ 1 03/20-3/2F(x), (2.15) 
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where we have defined the auxiliary functions 

y{x) := ^° T 2 (r) (x) , := 1 + y(x) . (2.16) 

^1^-103/20-3/2 

Since y possesses zeros on the real axis due to Uj and 9 h , we also define a renormalized 
function y^ 

M u N h 

y(x) = J]tanh^(x- Mj ) JJ tanh - h ) , (2.17) 

j=i h=i 

which obeys the functional relation 

y (r) (x - z'/2) y (r) (a; + i/2) = Q5 x (x + i'/2) ®i(ar - i'/2) , (2.18) 
as follows from fl2TT0l) and (1211]) . 



2.1.3 Derivation of NLIE 

The derivation of the NLIE can be most easily done in Fourier space. For a smooth function 
f(x), we define 



oo ^oo 



f[k] = — I e ikx f{x)dx, f(x)= I e~ ikx f[k}dk. (2.19) 



2vr 

We also introduce the special notation 

I r°° 

dlf[k] = — e ikx [\nf(x)}'dx (2.20) 

which will be frequently used below. 

It is convenient to introduce "shifted" Q functions, 

qi(x):=Q{x-i/2-ie), q 2 (x) := Q(x + i/2 + ie) . (2.21) 

By definition, qi is Analytic and NonZero (ANZ) for 3m x < 0, while q 2 is ANZ for Qm x > 0. 
We therefore have by Cauchy's theorem the important property 

dTq 2 [k >0] = dT qi [k<0] = 0. (2.22) 

Similarly, 

dhp a [k > 0] = dT(j) a [k > 0] = for a > , dhp a [k < 0] = d/0 a [A; < 0] = for a < . (2.23) 
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We slightly shift the arguments in (jglDj) . f[2TTTj) 

T 2 W (X + «C) = ^l+ e 01/ 2 +e03/2+ e ^""j ©(a?) 

g 2 (^ + i - lie) 



T^'\x-ie) = ^_i_ e 0_ 1/2 - e 0-3/2- 



Q3(x) . 



We then use the result 



1 

2^ 





hiT 2 (r) (x)" 







dx = i i 



(2.24) 
(2.25) 

(2.26) 



where we choose the contour C t as in Figure [fl and we obtain the following 



^ — 

X X— 






X 



le 



■le 



Figure 1: Integration contour 



— [fc > Q] = + + ^0-1/2-e + t«»[fc] - e" 2fce t«» [fc] 



1 + e~ k 

rdkdi. —ke 

Eg n 
1 + e~ k ' 



1 + e- 



(2.27) 



— r , nl <^i+ e + dlcj> 3/2+e + df0 1/2+e dZQ3[fc] - e 2k *dm[k] 
alq 2 [k < UJ — ; — ; — z r 



l + e k 



l + e k 



ikd h +kt 



+ e k 



(2.28) 



In addition, from (12.151) . one derives 

■fir \k\ - ^Tl + ^3/2 

U*! - e k/2 + g-fe/2 + 

where — (A; > 0) and < 0). 



jfc/2 _|_ e -fc/2 



e ifci? a +A;/2 
gfc/2 _j_ g-fc/2 



(2.29) 



8 



We shift the arguments in (12. 9p 

K x) = *&±fl_ r^-^ + fc), 

^1+603/2+601/2+6 - « + 2ze) 

sw - ^^.^^r-'L) ^ + > /2 - *> ■ « 2 - 3o > 

and then take the Fourier transformation. The substitution of (12.271 ), (12.281 ) and (12.291) 
into the resultant transformation then leads to the NLIE in Fourier space, 

ike h +ek e ik0 a +ek 

dib[k>o] = 1 + e _ fc +»LTT7 + i L e ^ + e -^ 



/i 



d/0 1/2+e , . ^ e ifcflh+£fc , ^ e ^+6fc 

dZ6[fc<o] = - -j^r- + » 2^ + * 2^ e k/2 + e ^ fc/2 

fc/2+efc 1 

+ e k/2 + e - fc/2 ^W + ^t^M ~ ^dm\k]) • (2-32) 

Interestingly, although a contribution from the inhomogeneities (ip) appeared during the 
calculation, it cancelled in the final form. An equation for y is immediately derived from 

(EH}, 

dly[k] = iJ^Y—^ + iJ^Y^k 



_ _|_ e « i _|_ c ft 

,fc(l/2-6) e -ft(l/2-e) 



+ efc /2 + e - fc /2 ^W + efc /2 +e - fc /2 ^W- (2-33) 

In the original coordinate space, the resultant equations read 

/oo /*oo 
G s (x - x') In «B(x') da/ - / G s (s - a/ + 2ie) In <8(a/) da/ 
OO ./ — oo 

/oo 
if(x-a/- z/2 + ie) In F(a/) da/, (2.34) 
-oo 

/oo 
if (a; -x' + i/2- ie) m <B(a/) da/ 
-oo 

/oo 
if (x - a/ - z/2 + ie) In «B(a/) da/ , (2.35) 
-oo 

where 

i /"oo — ikx i roo i i 

GJx) = — / rrrdk , if (ar) = — / ttttt^ dk = — — . (2.36) 

sK ' 2vr 1 + el fc l ' v; 2tt 7_ 00 2 cosh(A;/2) 2cosh(7rx) v ; 
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The source term in (I2.34p consists of the bulk ( "driving" ) contribution and the contribution 
from the hole excitations, 



D b (x) = Z^i(x) + £>£ le (z) - \N h , (2.37) 



where 



D$hW = L Xk(x) , £>2.(x) = Xk(x -$ a ) + J2 ^ - h) , (2-38) 

a=l h=l 

and 

X ' K {x) = 2ttK(x) , x'(x) = 2txG s (x) . (2.39) 

In particular, on suitable domains (containing the positive real axis), 

1 sinh(7r(x + i/2)/2) n n 

Xk(x) = -mtanh(7r(x - i 2) 2) = zln ... . . + — = arctan(smh(7rx)) - - , 

% smh(7r(x — 1/2)/2) 2 2 

(2.40) 

where Xk(Q) = —' n ~/2, and also x(0) = 0. The source term in (I2.35P is given by 

D y {x) = J2xK^-0h + i/2) + J2xK(x-Uj + i/2). (2.41) 

h=l j=l 

The parameters (uj, $ a , Oh) must actually be determined again by NLIEs. Indeed, (12.101) 
implies that the hole rapidities h are determined by 

b(9 h -ie) = bi(6 h ) = -l, (2.42) 

which also leads to the determination of the spinon-spinon and spinon-magnon scattering 
matrices, as discussed in Sec. HI 

In order to fix the parameters we need another NLIE. We consider the most natural 
auxiliary function a(x), defined bye 

X^x + i/2) X 2 (x-i/2) 

&\X) := "T — z TT = "J — 7 77^1 (2.43) 

V ; X 2 (x + i/2) X 3 (x-i/2) V ; 

where again A«(x) are defined in (12.51) . From (12 .4p we have 

Tl (x) = Ml 9t^A [i + a{x)] . (2.44) 



5 As discussed further in Sec. 12. 3[ one can verify numerically that ima(a;) and also 3? e [ilnb(x)] are 
increasing functions of x. 
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Hence, the zeros of Ti on the real axis $ a satisfy 

a(0 a ) = -l. (2.45) 

We omit the derivation of the NLIE for a(x) for |Qma;| < 1/2, which is similar to the 
one for the trigonometric and homogeneous case considered in [21]. The result is 

/oo roc 
K(x - x' - ie) In ©(a/) da/ - / K(x - x' + ie) In «B(x') dx' (2.46) 
-oo •/ — oo 

where the source term is given by 

N h M u 
h=l j=l 



2.2 The second set of BAEs ( Ej ) 



We finally consider an equation to fix the magnon rapidities Uj. For this purpose, we propose 
an expression for the transfer matrix eigenvalues similar to the one for the su(3) spin chain, 

/ \ -/ .sQ(x + i) ip(x + i/2) Q(x - i) .if>(x-3i/2) 
Tlx) = d> Ix — i) — — — V(plx + i)— — — — \-6{x — 2i)—- -— 

:= n{x) +t 2 (x) +t 3 {x). (2.48) 

Indeed, demanding analyticity of r(x) at x — lj (zeros of Q(x)) gives the BAEs ( 12.11) . while 
demanding analyticity at x = Uj + i/2 (zeros of i/;(x — i/2)) gives the BAEs ( 12. 2ft . 

Because of its similarity to the su(3) transfer matrix eigenvalue, we shall assume that 
t(x) is ANZ in the strip —1/2 < Qmx < 1/2, which is indeed the analyticity property for 
the su(3) case. This assumption can in principle be checked numerically for small values of 
L. However, we have so far not succeeded to do so, due to the difficulty of finding numerical 
solutions of the BAEs (El]), flZ2D - 

This assumption leads to a simple determination of Uj as follows. Let us consider an 
auxiliary function introduced in studies of the supersymmetric tJ model [2H] and the su(3) 
vertex model [29] . 

V ; n(x + i/2) +t 2 (x + i/2) v ; 



6 We expect that, starting from a suitable su(2, 1) R matrix, a transfer matrix can be constructed with 
eigenvalues (|2.48jl . However, we have not attempted to carry out this construction. 
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It is easy to check that this can be rewritten in terms of Ti(x) in (12.41) . 

<x) = « x -M : -*/2) . ( , 50) 



We then have 



which follows from 



c( Uj ) = -1 , (2.51) 



tlx) ~ 1 + c(x) - r(s + t/2) _ rji+jgjHj 

From our above assumption on the analyticity of t(x), the zeros of £(x) near the real axis 
are determined by those of ip (x) , namely Uj . 

The NLIE for c is obtained from the knowledge of T\. The result is 



/oo 
K(x-x' + i/2)\nY(x') dx' 
-oo 



(2.53) 



where the source term is given by 

M„ No 

D c (x) = Lx 2 {x) + ^2x3/2(x-u j )+^XK(x-'d a )-^(L + M u ), (2.54) 



2 

j=l a=l 



and 

1 f°° p—a-\k\—ikx 

X ' a (x) = 27rK a (x), K a (x) = —j -—^dk, Xa (0) = , a = 3/2 , 2 . (2.55) 
2.3 Counting functions and counting equations 

So-called counting equations relating the various types of Bethe roots and excitations in 
a given state can be derived from corresponding counting functions associated with the 
auxiliary functions . These counting equations help determine the spins of the excitations. 

We continue to restrict to the case of real w-roots and no r-roots. As in previous studies 
[211 [241 125]. it is convenient to classify Z-roots according to their imaginary parts as follows: 



2-strings : pairs of complex-conjugate roots Xj±iyj with < yj — 1/2 « 1, j = 1,..., M 2 /2 
real roots : Qrnlj — 0, j — 1, . . . , M rea i 
inner roots : \Qmlj\ < 1/2, j = 1, . . . , Mj 
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close roots : 1/2 < \Qm lj\ < 3/2, j = 1, . . . , M c 
wide roots : \^smlj\ > 3/2, j = 1, . . . , Mw 

Hence, 

M x = M real + M 2 + M I + M C + M W . (2.56) 

It is also convenient to introduce the functions 

„ , s ]- , f x — ia\ , 
T (x, a) = - In =f . (2.57) 



l \ x + xa 

Note that 9-(x, a) = 2 arctan(x/a) has branch points in the complex x plane at x — ±ia; 
following [21], we choose the corresponding branch cuts to be parallel to the real axis, 
extending from iat to +oo +ia, and from — oo — iat to — ia. This function has a discontinuity 
of — 2it when crossing the cuts from below. Similarly, we add to + (x, a) a 27r-discontinuity 
at x = so that it is a continuous function of x. 

We define the counting function z a (x) associated with the auxiliary function a(x) (12.431) 

by 

Mi M u 

z a {x) = -Loga(x) = L0_(x,l) -^^(x-^-,1) + s ^e + (x-u j , 1/2) . (2.58) 

We have verified numerically for various states that z a (x) is a continuous increasing function 
of x. This function "counts" zeros of T\(x) and real /-roots. That is, 

z a (xj) = 2TrI?, (2.59) 

where Jj 1 is integer (Si — S 2 odd) or half-odd integer (Si — S 2 even) if Xj is a zero of Ti(x) 
or a real /-root. Defining integers or half-odd integers I^ ax and I^ in by 

2; a (+oo) = 2tt (i^ + , 

z a (-oo) = 2tt - ij , (2.60) 

it follows from (12.581) and (I2.59j) . respectively, that 

Imax ~ Imin + 1 = Si + S 2 + M b 

= N» + M real , (2.61) 
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where is the number of /-roots lj with \Qmlj\ > 1. We therefore arrive at the first 
counting equation 

N 6 = S 1 + S 2 + M b - M real . (2.62) 
Similarly, we define the counting function Zf,(x) associated with the auxiliary function 

bt(x) by 

1 i r i i Mu 

zJx) = 9fte -Log bi(x) = 9fte \ - In lH — 7 — - + > 0+(a; — tt,-, 1) 

% Lz a(x — y2)J ^— ' 

Mi 



5=1 



+ L [0_(ar, 1/2) + 0_(z, 3/2)] - [0_(x - Z i5 1/2) + 0_(a; - Z i; 3/2)] } . (2.63) 

The presence of the first term generally requires the introduction of further discontinuities. 
We have verified numerically that Zb(x) is also a continuous increasing function of x. This 
function "counts" zeros of t!(\x) and centers of 2-strings and inner pairs. Proceeding as 
before, we find 



.„„,, I b mm + l = 2S 1 + M c + 2M w + (M 2 + M I )/2 

= N h + {M 2 + Mj)/2. (2.64) 
We therefore arrive at the second counting equation 

N h = 25i + M C + 2M W . (2.65) 

Finally, we define the counting function z c (x) associated with the auxiliary function c(x) 

(12321) by 

11 Mu 
z c (x) = Ke-Logc(x) =3tej - - In [1 + a(x + ?/ 2 )] + ^2^+(x - uj,l) 

M| 

+ L0_(x, 3/2) -J2 9 -( x - / i' 1 / 2 )}- (2.66) 

i=i 

We have verified numerically (using for the first term the same discontinuities introduced for 
the first term in (I2.63P ) that z c (x) is a continuous increasing function of x. Assuming 

IL* -ILn + ^ = L + M u - M real - (M 2 + Mj)/2 , (2.67) 

which can also be verified numerically, we recover the result 

M u = 2S 2 . (2.68) 
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3 Spin, energy and momentum of excitations 



We now compute the excitations' spin, energy and momentum, which enter into the compu- 
tation of the S matrix. Our results agree (except for some minor discrepancies) with those 
obtained previously using the string hypothesis. 

We can infer the spins of the excitations with the help of the counting equations found in 
Sec. 12.31 The second counting equation (12.651) implies that a spinon has Si = 1/2. Indeed, 
Nh = 1 requires S\ = 1/2 (and Mc = Mw — 0); Nh — 2 requires either Si — or Si = 1, 
etc. Note that all the terms on the RHS of (12.651) are nonnegative. Evidently, a spinon also 
has 5*2 = 0. The fact that a spinon has spin-1/2 was found using the string hypothesis by 
Takhtajan [llj. 

Similarly, the third counting equation f!2.68|) implies that a magnon has S2 = 1/2, and 
evidently Si = 0. This result was found using the string hypothesis by Beisert et al. [T5] . 

The spin of the $ particle is not determined by the first counting equation (12.651) . since 
not all the terms on the RHS are nonnegative. Nevertheless, an analysis of various examples 
suggests that this particle has Si = S2 = 0. 

By the definition in [13], the energy (E) is related to the anomalous dimension (jl.9p by 
7 = ^y^-E, and is therefore given by[] 

J= l 3 3 =l 3 

We can relate this to the derivate of the eigenvalue T 2 (x) (12. 5p at x = i/2, 



E = i—\nT 2 (x) 
dx 



3L 



M u 



OJ-J s. — ^ 



,=,/2 _ . =1 



/ I " 



_w| + 9/4 \uj-ij2 Uj-3i/2j_ 
Recalling the definition of the auxiliary function y(x) (12.161) . we see that 



(3.2) 



d Mu 



dx 

We observe from ( 12.20p that 



E = i—\ny(x) - 2L + V ( ; = -J. (3.3) 



/oo 
dk e~ ikx dly[k] , (3.4) 
-00 



d_ 

dx 



7 For convenience, we drop the constant term 7L/6 in the expression for E. This definition of energy is 
(for the Z-roots) a factor 2 larger than the one in pTj . 
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and substitute our result for dly[k] f!2.33j) to obtain 

3 1 



E 



M u 

3=1 



7T 



+ 



cosh(7TM j ) u) + 9/4 u) + 1/4 



^ COSh(7T0fc) 



where the ellipsis (...) represents the Casimir energy contribution. We conclude that the 
energy of a spinon is 



e h {9) 



7T 



cosh(7r6 ) ) 



and the energy of a magnon is 



7T 



+ 



(3.6) 



(3.7) 



v ' COsh(TTM) M 2 + 9/4 M 2 + l/4' 

in agreement with Eqs. (6.15), (6.32) in Beisert et al. [13] . respectively, up to a factor 2. 
The spinon result (13.61) was first found by Takhtajan [IT] . We remark that 



e u (tt) = 2ttK 2 {u) , 



(3i 



where K2(u) is the kernel introduced in (12.551) . Evidently there is no ^-dependent contribu- 
tion in (I3.5p . which implies that the $ "particle" does not carry energy. 

The momentum is given by % 

Mi /, , -\ M u 



p = i 



5> 

.3=1 



+ £ln 



uj — 3i/2 
itj + 3i/2 



(mod 27r) 



(3.9) 



We can evaluate it in similar fashion. Indeed, we find that 

Mu 



P = -lny(i/2) + - ^ [lne_ 3 (%) + lnei(i^)] + L7r 



3=1 



where we have introduced the notation 



e n \u 



u + in/2 
u — in/2 

Proceeding as before, we arrive at the result 

M u 

P = + [Xk(Uj) + q 3 (uj) - ft («,■)] + ^Xk{0h) + ■■■ , 

3=1 h 

where Xk( x ) is defined in (12.391) . and q n (x) is defined by 

q n (x) = n + i\iae n (x) n>0, q_ n (x) = -q n (x) , q (x) = 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



3 This definition of momentum differs (for the i-roots) by an overall sign from the one in 
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It is an odd function of x, and satisfies 



Qn( x ) — 2 arctan (2x/n) n^O. (3-14) 
We conclude that the momentum of a spinon is 

Ph{0) = XK(0), (3.15) 

and the momentum of a magnon is 

Pu(u) = Xk{u) + q 3 (u) - qi(u) , (3.16) 

in agreement with Eqs. (6.15), (6.32) in [13], respectively, up to an overall sign. Correspond- 
ing to the energy result (13.81) . we observe that 

Pu(u)=X2(u), (3.17) 

where X2{u) is defined in (12.551) . The particle also does not carry momentum. 

4 S matrix 

We finally turn to the problem of computing the scattering amplitudes for the various exci- 
tations. 

4.1 Spinon-spinon 

It is convenient to review the computation of the spinon-spinon S matrix [181 OS] using the 
NLIE approach [25]. Let 9^ , 9h 2 denote the rapidities of the two spinons. Since b(9f ll — ie) = 
-1 (T222D, the lnb equation f[2T34j) implies 

/oo 
dx' K(9 hl - x' - i/2) In Y(x') , (4.1) 
-oo 

since the convolution terms involving 03 and 23 become exponentially small in the IR limit. 
Neglecting the convolution term in the Iny equation (I2.35p . one obtains 

y(x) = tanh(7r(x-^ 1 )/2)tanh(7r(x-^ 2 )/2), (4.2) 

and therefore 

Y(x) = 1 + v(x) = cosh(n(x-(9 hl +9 h2 )/2)) 

{ ) i + cosh(7r(x-^ 1 )/2)cosh(7r(x-^ 2 )/2) ' { ' 
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We now exponentiate both sides of (14.11) . and note using (j2.37p . (12.381) that 



D b (6 hl ) = L X K(0 hl ) + x(0) -K, e = 6 hl -9 h2 . (4.4) 

With the help of the momentum expression (13.151) . we compare the result with the Yang 
equation 

e iL Ph (e hl ) s Kh (0) = l. (4.5) 
We conclude that the 5* matrix is given (up to a constant) by 

S h M = e i xWS RSOS (e), (4.6) 



where 



C (a\ ^f°° dx' K(9 h ,-x'-i/2)lnY(x') 

= e -^Ww(«)-vaW] _ e -iiM8)-M0)\ ? (4.7) 



and 



sinh(7r(« + x) A) 

ipo(x) = arctansinh(7ra:/2) = im . , , tttt , (4.8) 

smh(7r(2 — x) /4) 

_ f°°„ sjn(fcx)sinh((n-l)fc/2) 
<p n {x) - ^ UK A , sinh(nA;/2)cosh(A . /2) > l 4 - y J 

with v?4(x) = (<^2(^) + V'o^)) /2- The convolution integrals are performed using the results 
collected in the appendix. The result ( 14. 7ft is (up to a crossing factor, and a rescaling of 
the rapidity by n) one of the kink-kink scattering amplitudes of the tricritical Ising model 
perturbed by the operator $(1,3) [IT] , which appears also in the soliton-soliton S matrix of 
the supersymmetric sine-Gordon model [IB] . We note that 

1 r(i + »fl/2)r(i/2-tfl/2) 
xW -7 ln r(i-^/2)r(i/2 + ^/2)' (4 - 10) 

which is (up to the same rescaling of the rapidity by 7r) the soliton-soliton scattering phase 
of the sine-Gordon model [30] in the isotropic limit j3 2 — > 87r. 

If we also consider the In y equation with an additional Z7r term, then the RHS of (14.21) 
acquires a minus sign. The corresponding amplitudes can be computed along similar lines 
|25j . However, for simplicity, we restrict our attention to the \ay equation without this 
additional m term. 
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4.2 Spinon-magnon 

Let 9^ , «i denote the rapidities of the spinon and magnon, respectively. The spinon-magnon 
S matrix can be computed in two different ways. One way is to start from b{9y Ll —ie) = — 1, 
which again leads to (14. ip . The lay equation implies 

y(x) = tanh(7r(x - 9 hl )/2) tanh(vr(x - u%)/2) , (4.11) 

and therefore 

Y(x) = 1 + y(x) = cosh(rr( X - (9 hl + Ml )/2)) 

K ' yv ' cosh(7r(a:-0 ftl )/2)cosh(7r(z-ui)/2) v ; 

Moreover, now D&^J = LxK(@hi) — -^Pa(^Ai)) u P to an additive constant. Proceeding as 
before, we obtain the result 

Sm(0) = ^«so5(^) , (4.13) 

where now 9 = 9 hl — Ui, and Srsos(9) is given by (14. 7p . That is, in contrast to the spinon- 
spinon 5 matrix (14. 6ft . the spinon-magnon 5 matrix consists only of the RSOS factor. 

A second way to compute the spinon-magnon S matrix is to start from c(«i) = —1 (I2.5ip . 
which together with the In c equation (I2.53P imply 



m = iD c ( Ul ) + / dx' K{u x -x' - i/2)lnY(x') . (4.14) 

J — oo 

We exponentiate both sides of this equation, and note that 

A,(ui) = Lx2{ui) = Lp u (ui) , (4.15) 

where we have made use of (I2.54p and the momentum result (I3.17p . Comparing with the 
corresponding Yang equation, we recover the same result, i.e. 

S u>h (9) = S RSO s(0) , (4.16) 

where now 9 = u\ — 9y ll ■ 
4.3 Magnon- magnon 

Let Mi, ui be the rapidities of the two magnons. The lny equation (12.351) implies 

y(x) = tanh(7r(x — Mi)/2) tanh(7r(x — u 2 )/2) , (4-17) 
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and 

v ; yy ' cosh(7r(2;-iti)/2)cosh(7r(a;-w 2 )/2) v ; 

The condition c{u\) = — 1 (I2.5ip and the lnc equation (12.531) again give (14.141) . where now 

(cf. dnsj) 

D c (u 1 ) = Lp u (u 1 )+Xa/2{0), (4.19) 

with 9 = u\ — U2- Proceeding as before, we conclude that the magnon-magnon S matrix is 
given by 

S u ,u{9) = e lx ^ {e) Srso S (6) , (4.20) 

where Srsos(9) is given by (14.71) . We note that 

1, T(-l/2 + i9/2)T(-i9/2) 
v 3/2 (fl) = -ln \ ' — 1 ,\ \ — f-^ + vr, 4.21 

and that s(9) = e iX3 / 2(e) has the crossing property 

'\-%9 
2 + i9 

Hence, s(9)/(l +i9/2) is crossing invariant. 

We have considered so far the composite operators containing only D a \ covariant deriva- 
tives and computed the S matrix amplitude between them. In principle, one would need to 
add r-roots to compute amplitudes for the derivatives carrying the right-spin state 2. But 
this can be done, without adding r-roots, by using the SU{2) R symmetry. The "vertex" 
part of the S matrix is in fact a 4 x 4 matrix which can be fixed completely by the SU(2)r 
symmetry along with factorizability (i.e., Yang-Baxter equation), unitarity and crossing, 

s(9) 



s (i _ 0) = ( LJ? ) s (0) . (4.22) 



1 + i9/2 

where V is the permutation matrix. 



(V + i0/2) , (4.23) 



4.4 ^-spinon and ^-magnon 

The condition fl(i? a ) = — 1 (I2.45P together with the lno equation (I2.46P imply that the S 
matrices S$ t h and S# )U are identical, and are given by 

_ siDh(7r((9/2-i/4)) 

m ~ ^MMm+m ' ( } 

The same result can also be obtained starting from (12.421) . (I2.34p (for Sh,&) and from (12. 5 1 j) . 
(I2.53P (for S u> 0). Since there is no ^-dependent contribution in the source term of the In a 
equation (12.461) . there is no nontrivial •d-'d scattering. 
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5 Discussion 



We have proposed a set of NLIEs (12331) - (12^11) . (I2T46D . (1233), (EJ^-tES!) to describe the 
QCD spin chain of Beisert et al. [13]. We have used these NLIEs to compute S matrix 
elements for excitations of this model, as shown in detail in Sec. HI The consistency of our 
results (S a ,b = Sb, a for particles a and b of different types) provides further support for the 
validity of these NLIEs. 

Many questions remain to be addressed. It should be possible to generalize this work 
along the lines [31] and compute the boundary 5* matrix for the open QCD spin chain 
corresponding to operators with quarks at the ends. The magnon-magnon S matrix (I4.20p . 
(I4.2ip has an infinite number of singularities (starting at 9 = ±2i), which can presumably 
be interpreted as magnon-magnon bound states ("breathers"). The energy and momentum 
of these breathers was computed using the string hypothesis in [13]. It would be interesting 
to analyze these excitations without invoking the string hypothesis, and to determine their 
S matrices. It would also be interesting to consider the effects of higher loops ([7] and 
[T3] worked only to leading order in the 't Hoof coupling) and to better understand the 
significance of these results for QCD, as well as for the full M = 4 SYM theory and for the 
corresponding string theory. 



Acknowledgments 

One of us (CA) thanks Shizuoka University and the University of Miami for support. This 
work was supported in part by a Korea Research Foundation Grant funded by the Korean 
government (MOEHRD) (KRF-2006-312-C00096) (CA), by the National Science Foundation 
under Grants PHY-0244261 and PHY-0554821 (RN), and by the Ministry of Education of 
Japan, a Grant-in-Aid for Scientific Research 17540354 (JS). 



21 



A Convolutions 

The convolution integrals involving the kernel K(x — |) = 2 S mh nx can ^ e evama t e d using 
the following results 

if 00 ,, lncosh(7r(x' — ie)) i . . . 1, , . . . 

- / ax — — — ; r- = — arctansinh(7ra;H — lncosh(7rx), (A.l) 

2J-00 smh(7r{x -x' + ie)) 2 v ; 2 v ; ' v ; 

1 f°° , , In cosh(7r(x' — ie)/2) % , . 1, , , . . .. . 

-/ dx' — - , ) v — ^-f = — (p 2 (x) +-lncosh(7rx/2), A.2 

2J-0O sinh(7r(x-x' + ie)) 2^ v ; 2 v 7 ; ' v ; 

if 00 -,, In sinhfvrfx' — ie)) 1, . , , , 1, _ . . s , . „ s 

-/ dx' — — JJ - = -lnsinh(vrx) - -lntanh(7rx/2) , (A.3) 

2 smh^x-x' + ie)) 2 v ' 2 v 7 7 4 ' v ' 

if 00 ,, lnsinh(7r(x' — ie)/2) i . . 1, , , . . . . 

-/ da/ — : . , ; — ^rv = (x + -lncosh vrx/2 , A.4) 

2 7_ 00 sinh(7r(x-x' + ie)) 2^ A ; 2 v 7 ; 4' v ; 

where e is a small positive number, and (p2(x) is given by ( 14. 9p . 
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